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Abstract 

The K2K experiment has presented the first results on the observation of u^. They 
show a depletion compared to the expectations and are consistent with neutrino oscilla- 
tions with a mass-splitting in the range favoured by the Super-Kamiokande atmospheric 
neutrino measurements. Here we examine the extent by which the range of Am^ obtained 
from the K2K measurements can vary due to the uncertainties in the flux, cross-section, 
and detector efficiency. 
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I Introduction 



A non-zero neutrino mass has far-reaching implications in particle physics, astrophysics, and 
cosmology. Any evidence in its support must be probed from as many different angles as 
possible. Atmospheric neutrinos provide a strong case for non-degenerate neutrino masses 
and mixing. Super-Kamiokande (SK) [1], MACRO [2], and Soudan 2 [3] studying neutrinos 
produced in interactions of cosmic rays with atmospheric nuclei show a significant suppresion 
of the ratio of to z/g with respect to the expectation from standard hadronic shower models. 
These observations can be explained by the phenomenon of neutrino oscillations in which 
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changes to some other neutrino, e.g., v^. Such oscillations can occur if neutrinos have a non- 
degenerate mass spectrum and if the neutrino mass eigenstates are not the same as the flavour 
eigenstates. The atmospheric neutrino data favour a mass splitting Am^ = — m|| ~ 10"^ 
eV^ and maximal mixing, ^ = 7r/4. Neutrino oscillations have also been advanced as a natural 
solution to the long-standing solar neutrino problem, especially when the data from the SK 
and SNO experiments are taken together. 

K2K [4] is the first accelerator based long-baseline oscillation experiment. It is designed pre- 
cisely to probe the mass and mixing parameters favoured by the atmospheric neutrino data. It 
uses neutrinos produced by the 12 GeV proton beam of the KEK accelerator laboratory which 
are subsequently detected by the SK detector at a distance of 250 km. The average neutrino 
energy is {E) ~ 1.3 GeV and the intense beam is 98.0% pure v^^. K2K also uses a near detector 
at KEK, 300m downstream from the pion production target, to monitor and study the neutrino 
beam. The experiment started taking data from June, 1999. The primary oscillation search 
mode is v^, —>■ oscillation (the disappearance mode). 

The first results from K2K [5] have been released. With 3.85 xlO^^ protons on target (p.o.t.), 
in the 22.5 kton fiducial volume of the SK detector they observe 44 fully contained events in 
place of the (no oscillation) expectation of 63.9lg;e events. This corresponds to a depletion 
ratio 

^ = No. of observed events ^ ^ ggg+o.OTQ. 
No. of expected events -o.oeo 

This depletion is consistent with neutrino oscillations with a mass splitting in the atmospheric 
range and maximal mixing^. In this note, basing ourselves on the available information, we 
wish to examine to what extent this conclusion could be uncertain. Our finding is that the 
result is robust. 

When oscillations are operative, the survival probability of a neutrino with energy E produced 
a distance L from a detector is expressed (in the two-flavour approximation) as: 

^ ^, X 2.2 1.27Am2(eV^)L(km) 

Pmm = P^^^ ^ ^m) = 1 - sin' 2^ sin^ ^^^^^^j ^ ^ (2) 

For the purpose of the present work, the nature of the partner neutrino state to which the v^J^ 
oscillates is immaterial so long as it does not mimic a at the SK detector. 



II The K2K results 

In the presence of neutrino oscillations, the event rate in the K2K experiment, r, can be 
expressed as 



^Based on 2.6x10^^ p.o.t. [6] the depletion ratio was d = 0.670lg:2^^- 
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r = I dES{E)P,,{E), (3) 

where P^^ is given by (2). 

The oscillation parameters 9 and Am^ will be determined from the K2K data. The experiment is 
sensitive to (Am^,sin^2^) around (10~^eV^, 1) favoured by the SuperKamiokande experiment. 

S{E) is the 'energy spectrum' of the neutrinos from KEK as detected at the Kamioka site. This 

spectrum is given by the convolution of the standard (no oscillation) neutrino energy spectrum 
at the SK (far) detector, the interaction cross-section a, and the detection efficiency e (all in 
differential form), integrated over the final state parameters, i.e., symbolically, S = J (pae. 

In the published literature there is information on (j) but not yet on the product a.e relevant 
for the far detector specification. As direct and accurate reconstruction of the spectrum is not 
possible, some indirect, approximate reconstruction was suggested by the K2K collaboration. 

We have two sets of information in hand. 

(i) The neutrino spectrum is considered to be practically zero above 5 GeV [4]. 

(ii) According to the K2K Monte Carlo simulation, for an exposure of 3.85 x 10^^ p.o.t. , assuming 
maximal mixing, we have the set of predictions for the event rate (r^*^) for different Am^ shown 
in Table 1: 



Am^ in eV^ 


1 X 10-^ 


3 X 10~3 


5 X 10-3 


7 X 10-3 




63.9 


41.5 


27.4 


23.1 



Table 1: The values of r^'^ for different Am^ determined from the K2K Monte Carlo calculations [5] for 
maximal mixing (sin^ 20 = 1). 



Ill Impact of different *S'(^) 

In the absence of direct information, a form for the spectrum was suggested [7] to be S{E) oc 
x°(l — xy with X — E/{5 GeV). It fits the expected behaviour of S{E) at x — and 1 and 
has a particularly simple form. The overall normalization and the parameters a and f3 were 
obtained by using a best fit to the data in Table 1. It was shown using this functional form 
of S{E) that the preliminary data from K2K are consistent with the Am^ favoured by the SK 
atmospheric neutrino experiment. 

Here, we ask the question: How sensitively does this conclusion depend on the form of S{E)'! 
To this end, we consider several different functional forms for the spectrum S{E). For each, 
the best-fit values of the parameters are determined by minimizing 

= E ^ ' J ^ , (4) 

i=2 '^i 
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where is calculated using the chosen form of S{E) and rf^*^ is taken from the K2K Monte 
Carlo results of Table 1 and, as in [7], 

= 44 + O.Olr^ (5) 

r*^*" corresponding to i = 1, i.e., Am^ = 1.0 x 10~^ eV^, is used to determine the overall nor- 
malization constant N. In this way, in all cases the no oscillation limit is correctly reproduced. 
The functional forms that we have examined and the best-fit values of the parameters are listed 
in Table 2. These functions are plotted in Fig. 1. 



Case 


functional 
form 


Normalization 

N 


best fit values 


a 


/9 




a 




1370.7 


1.14 


2.78 


0.044 


b 




518.3 


0.85 


3.81 


0.063 


c 


Nx'^ll - x^y 


307.7 


0.60 


5.50 


0.068 


d 




215.7 


0.42 


7.48 


0.061 



Table 2: The different parametrizations of S{E), the best-fit values of the parameters, and the values of 

x'. 

All the above functions have a common nature in that they vanish at a; = and 1 and rise to a 
maximum value at some intermediate energy. The actual neutrino spectrum as a function of the 
energy is expected to have rising and falling regions at the two ends and an intermediate region 
where its variation is rather slow. To mimic this feature, we have also considered two further 
cases where we have chosen the functional forms of Cases a and b of Table 2 but 'chopped' the 
function at the values of E where it achieves half the maximum value. Between these points. 
El and E2, the function is assumed to be a constant. The points Xi — Ei/ (5 GeV), (i =1,2) and 
the best-fit values of the parameters in these cases are presented in Table 3 and the functions 
exhibited in Fig. 1. 



Case 


functional 
form 


Normalization 

N 


Half-Maxima 


best fit values 


Xi 


X2 


a 




x' 


e 


Nx'^il-xf 


13136.3 


0.128 


0.529 


2.00 


4.61 


0.065 


/ 


Nx^'il-xY 


3392.6 


0.156 


0.544 


1.72 


6.58 


0.047 



Table 3: Parametrizations of S{E) with a constant region between xi and X2 (see text). The best-fit 
values of the parameters and the values of are shown. 

We expect that the six rather different functional forms for S{E) that we have chosen adequately 
capture the possible uncertainties in the neutrino spectrum, cross-section, and detection effi- 
ciency. We make two remarks. Firstly, the fits are all of very low x^- This is due to the 
comparatively large af from (5) associated with each datum point. Secondly, we find that over 



4 



a rather broad region in the parameter space around the best-fit points, varies rela- 

tively little. In particular, the best-fit point we find in Case a is somewhat different from the 
one in [7] even though the functional forms are the same. For the latter we find a comparable 
but higher x^. 

After having obtained the six different functional forms of S{E) which best fit the Monte Carlo 
results of Table 1 we use them in (3) to obtain the predicted value for the number of events for 
the K2K experiment with 3.85 x 10^^ p.o.t. These results are shown in Fig. 2 as a function of 
Am^ for the maximal mixing case (sin^ 26 = 1). The two horizontal lines are the la range of the 
experimental observation. The intercepts of the curves with this experimental range determine 
the allowed values of Am^ for the various cases. It is seen from Fig. 2 that the differences 
between the alternative cases (a - /) are comparatively marginal. There is a small spread in 
the lowest admissible values of Am^ in the different cases. But as the data improve even this 
will be removed. Since we have chosen S{E) with a wide range of forms, it would not be unfair 
to conclude that the range of Am^ which is favoured by the K2K data can be predicted in a 
robust fashion. 

Motivated by the best-fit to the atmospheric neutrino data, in this work we have used maximal 
mixing, i.e., siv? 29 = 1.0, in (2). The SK data restrict the mixing angle rather tightly and at 
90% C.L. one can allow sin^ 20^0.9. How much would our conclusions be affected if the lower 
limit of sin^ 29 is used in the analysis? We find, for example, for an S{E) of the form of Case 
a (see Table 2) the best fit point moves to {a,f3) = (2.39, 6.10)^. The event rate as a function 
of Am^ for this case is also shown in Fig. 2. Notice that a small region around 2 x 10~^ eV^ is 
allowed for this case. A large improvement in the K2K systematic and statistical uncertainties 
will be able to exclude this option. 

IV Conclusion 

The first results from the K2K long baseline experiment show a depletion in the number of 
events compared to the expectation. This is consistent with oscillations of the z/^ to a neutrino 
of a different fiavour with mass and mixing in the range favoured by the atmospheric neutrino 
results. In this note, we have shown that this conclusion is robust and is not altered even if the 
uncertainty in the initial neutrino spectrum, interaction cross-section, and detection efficiency 
are widely varied. 

Note added: After the submission of the paper for publication, K2K announced results [8] based 
on 4.8 X 10^^ p.o.t. They observe 56 events with a no-oscillation expectation of 80.ll5;4. This 

corresponds to the depletion ratio (see Eq. (1)) d = 0.699lo;o5o, consistent with the earlier result. 
We find using our analysis that for the different paramctrizations of S{E) this corresponds to 
allowed ranges for Am^ shown in Table 4. It is gratifying that the best fit point found by the 
K2K group for maximal mixing (Am^ = 2.8 x 10~^ eV^ [8]) lies within the allowed range for 
all parametrizations considered. 

^This fit has = 9.74 x IQ-^. 
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Case 


Am^ Limits (10"^ eV^) 




Lower 


Upper 


a 


1 ^7 


fil 




1.(30 


3.72 


c 


1.52 


3.40 




22.48 


25.06 


d 


1.47 


3.73 




21.77 


26.38 



Table 4: The ranges of Am^ allowed at la for maximal mixing by the latest K2K results[8] for the 
alternative parametrizations of S{E) (see Table 2). Notice that in cases (c) and (d) there are two 
allowed ranges. 
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Figure 1: The spectrum S(E) as a function of E. The cases to which the different curves correspond 
are indicated (see Tables 2 and 3). 
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Figure 2: The number of events in the K2K experiment for an exposure of 3.85 x 10^^ p.o.t. as a 
function of Am^ (at sin^ 2^ = 1). The region between the two horizontal hues is the current K2K 
measured range within 1 standard deviation. The boxes are the values from the K2K MC simulation 
(Table 1). The six curves obtained from the different forms for S{E) considered in this work arc shown 
(same conventions as in Fig. 1). The curve marked 'non-max' corresponds to the case of non-maximal 
mixing (sin^ 29 = 0.9) discussed in the text. 
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